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DYNAMICS  OF  LIQUID-FILLED  SHELL: 

RESONANCE  AND  EFFECT  OF  VISCOSITY 

ABSTRACT 

Experiments  have  been  conducted  with  liquid-filled  shell  in  free 
flight  and  with  a  liquid-filled  gyroscope.  The  cavities  in  each  were 
designed  to  lave  a  principal  fluid  frequency  which  was  equal  to  the 
nutational  frequency  of  the  respective  systems. 

Comparison  of  experimental  results  with  predictions  of  the  inviscid 
theory  of  Stewartson  showed  significant  Reynolds  number  effects.  Both 
the  rate  of  divergence  of  the  nutational  amplitude  at  resonance  and  the 
width  of  the  resonance  band  are  Reynolds  number  dependent.  Also,  there 
is  a  snail  shift  of  resonance  frequency  with  Reynolds  number-  The 
experiments  suggest,  however,  that  in  the  limit,  for  large  Reynolds 
numbers  of  the  order  of  10^,  the  inviscid  theoiy  gives  excellent 
predictions . 
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1.  INTRODUCTION 


1* 

Earlier  experiments  in  the  free  flight  range  have  provided  data  on 
the  dynamic  "behavior  of  liquid-filled  shell  during  the  very  early  part 
of  their  trajectory.  Observation  showed  that  during  spin -up  of  the 
liquid,  the  shell  were  dynamically  unstable  for  all  fill  conditions  of 
practical  interest.  Moreover,  the  instability,  or  the  rate  of  divergence 
of  the  nutational  component  of  yaw,  depended  markedly  on  the  specific 
gravity  of  the  liquid;  i.e.,  for  heavier  liquids  the  divergence  of  yaw 
was  more  rapid.  For  example,  mercury-filled  shell  developed  more  than 
80°  of  yaw  in  150  feet  of  travel  from  the  muzzle.  It  was  also  found  that 
the  spin-up  of  the  liquid,  as  measured  by  the  loss  of  angular  momentum  of 
the  shell  in  vacuum,  was  markedly  accelerated  by  the  development  of 
secondary  flow  within  the  cavity  due  to  end  effects.  A  theory  for  this 

nhenomenon  has  been  developed  and  is  in  good  agreement  with  the  experi- 

.  2 
menus  . 

In  the  above  experiments,  no  specific  attempt  was  made  to  verify 
Stewartson’s  theory^  because  for  the  particular  cavity  used,  the 
predicted  resonance  occurred  only  at  45  percent  fill-ratio  and  relatively 
few  firings  were  done  at  this  loading. 

To  the  best  of  our  knowledge,  Stewartson's  theory  has  been  previously 

tested  only  twice.  The  first  experiments  were  done  by  firing,  at  long 

ranges,  specially  modified  shell  some  of  which  were  deliberately  designed 

to  be  unstable  (resonate)  according  to  Stewartson’s  criteria,  and  others 

to  be  stable.  These  experiments  were  spectacularly  successful.  The 

3 

second  experiments  were  made  by  G.  II.  Ward  with  a  gyrostat.  Ward 
found  good  agreement  with  the  theoretical  prediction  of  the  fill-ratio 
at  which  the  principal  mode  of  instability  should  occur.  However,  he 
also  found  that  the  resonance  band,  or  the  range  of  fill-ratio  over 
which  the  gyrostat  was  unstable,  was  more  than  twice  as  broad  as  theo¬ 
retically  predicted.  Considerations  of  various  sources  which  might  have 
contributed  to  this  discrepancy  led  to  inconclusive  results. 


Superscript  numbers  denote  references  found  on  page  39. 
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In  view  of  the  importance  which  Stewartson’ s  theory  nay  "be  expected 
to  play  in  the  design  of  liquid-filled  shell,  it  appeared  desirable  to 
examine  this  discrepancy  further.  Therefore,  free  flight  tests  were 
conducted  in  the  inclosed  range  under  carefully  controlled  experimental 
conditions.  These  tests  were  supplemented  by  tests  with  a  liquid-filled 
gyroscope.  The  gyroscope  permitted  tests  over  a  wide  range  of  Reynolds 
numbers . 

This  report  describes  the  results  of  these  experiments. 

2.  FREE  FLIGHT  FIRINGS 

The  shell  for  these  experiments  were  20mm  models.  A  schematic 
drawing  of  the  model  is  shown  in  Figure  1.  The  models  were  machined 
from  solid  aluminum  bar  stock  to  have  cylindrical  cavities  of  prescribed 
fineness  ratios.  The  cavities  were  designed  so  that  one  of  the  principal 
frequencies  of  the  fluid  was  close  to  the  nutational  frequency  of  the 
shell.  The  aerodynamic  characteristics  of  this  configuration  were 
well  known  from  previous  work  .  These  and  other  characteristics  as 
inferred  from  the  present  firings  are  discussed  in  Appendix  1. 

To  test  Stewartson’ s  theory,  the  fluid  must  be  fully  spinning. 
Consequently,  with  inclosed  range  firings  the  fluid  must  come  to  full 
spin  very  rapidly.  Previous  tests  showed  that  glycerine,  v  =  1000 
centistokes,  reaches  full  spin  in  about  40  feet  of  travel  from  the  muzzle. 
The  remaining  240  feet  of  the  range  could  then  be  used  for  relevant 
observations  during  which  the  shell  executes  about  ten  nutational 
oscillations.  From  prior  physical  and  aerodynamic  measurements  the  non- 
dimensional  nutational  frequency  of  the  model,  t  ,  was  estimated  to  be 

n  2  2 

0.06.  The  cavity  was  to  be  filled  to  90  percent  (b  /a  =  .10).  There¬ 
fore,  from  Stewartson’ s  tables^5'’,  any  cavity  with  height  2c  and 

diameter  2a  or  fineness  ratio 

c/a  =  1.051  (2j  +1)  3  =  0,  1,  2,.... 

will  have  as  its  principal  fluid  frequency  equal  to  the  nutational 
frequency  of  the  shell  =  .06.  For  3  =  1,  c/a  =  3*153* 
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In  order  to  cover  a  broader  range  of  fluid  frequencies,  it  was 
decided  to  set  the  fill-ratio  at  90  percent  and  vary  the  fineness  ratio 
from  3.03  to  3-29*  This  variation  of  c/a  permits  testing  of  shell 
response  to  fluid  frequencies  varying  from  t  q  =  .03  to  iQ  =  .09*  At  a 
constant  fineness  ratio,  say  c/a  =  3-153,  hy  varying  fill-ratio,  the 
lowest  attainable  frequency  at  100  percent  fill  would  have  been  to  =  .047 
producing  an  asymmetrical  distribution  about  t  n  =  -C6. 

The  observed  nutational  yaw  damping  rates  of  the  liquid-filled  models, 
corrected  for  normal  aerodynamic  damping,  are  shown  in  Figure  2.  These 
were  obtained  by  the  usual  analysis  of  range  firings.  Also  shown  are 
theoretically  predicted  rates.  As  can  be  seen,  there  are  significant 
differences  between  the  theoretically  predicted  and  the  observed  behaviors. 
We  shall  defer  discussion  of  these  differences  until  after  presentation 
of  the  experimental  results  from  the  gyroscope  tests. 

3.  GYROSCOPE 

3*1  Description  of  the  Apparatus: 

The  instrument  shown  in  Plate  1  was  originally  designed  and  built 
during  the  last  war  to  study  the  dynamics  of  liquid-filled  shell.  The 
rotor  was  a  modified  4"  shell  which  was  driven  by  an  air  turbine,  the 
air  being  admitted  through  the  outer  and  inner  gimbals. 

For  our  purpose  it  was  altered  by  discarding  certain  components. 

A  new  rotor  was  built  containing  a  removable  lucite  cylindrical  cavity 
of  maximum  inner  diameter  2.50"  and  height  12.90"  or  fineness  ratio 
c/a  =  5.16.  The  full  capacity  of  this  cavity  is  1038  cc.  Unfortunately 
the  center  of  the  cavity  was  above  the  pivot  point  which  presented 
certain  experimental  inconvenience.  However,  the  center  of  mass  of  the 
whole  system  could  be  varied  by  adjustable  rings  on  the  rotor  and  the 
inner  gimbal.  Most  of  the  experiments  were  conducted  with  the  gyroscope 
in  a  stable  position;  i.e.,  with  the  center  of  mass  of  the  system  below 
the  pivot  point.  The  inner  and  outer  gimbals  were  mounted  on  ball 
bearings . 
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The  oscillations,  about  one  gimbal  axis,  vere  measured  with  strain 
gages  mounted  on  a  strip  of  .015"  stainless  steel  shim  stock.  The  strip 
was  spring  loaded  and  attached  to  the  outer  gimbal,  see  Hate  1.  The 
gages  formed  part  of  a  bridge  circuit  whose  amplified  output  was  monitored 
on  a  photographic  recorder.  Two  records  are  shown  in  Figure  5- 

The  spin  of  the  rotor  was  measured  by  a  small  pick-up  coil  mounted 
on  the  inner  gimbal  and  two  small  magnets  on  the  rotor.  The  output  was 
recorded  on  a  tachometer  and  on  an  electronic  counter.  The  rotor  could 
be  spun  in  excess  of  10,000  rpm  with  reasonable  safety;  however,  most  of 
the  experiments  were  conducted  at  lower  spins. 


The  transverse  moment  of  inertia  of  the  gyroscope  could  be -varied, 
within  narrow’  limits,  by  the  adjustable  rings.  An  example  of  the 
characteristics  of  the  gyroscope  with  empty  cavity  (c/a  =  5-l6)  and 
certain  position  of  the  rings  is: 


Weight,  rotor  and  inner  gimbal,  lbs.  25-80 

2 

Moments  of  inertia:  ,  lbs -in  32.5?- 


...  -  2 
r  ,  Xbs-in 

865-5 

0 

.  2 

,  ,  lbs -on 

833-8 

0 

850 

=  .058 


3.2  Experimental  Results: 

From  preliminary  experiments  with  a  completely  filled  (c/a  =  5-l6) 

cavity,  the  nutational  frequency  of  the  gyroscope  was  found  to  be  closer 

to  .01*0  than  to  .038  as  given  above.  For  this  cavity,  from  Stewartson’s 

tables,  one  finds  that  resonance,  i.e.,  t  =  t  =  .(AO,  should  occur  at 
7  on 

92  percent  fill-ratio.  The  resonance,  at  this  fill-ratio,  should  occur 
in  any  cavity  whose  fineness  ratio  satisfies  the  condition 
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c/a  =  1.032  (2j  +1)  j  =  0,  1,  2.... 

Tims.,  'or  j  =  c/a  =  3*10;  for  j  =  2,  c/a  =  5-16,  etc.  Both  of  these 
cavities  were  used. 

In  order  to  cover  a  broad  range  of  Reynolds  lumbers,  the  following 
liquids  were  tested: 


Liquid 

0,  gm/’cc 

V,  c.s 

Mercury 

13.6 

0.1 

Water 

1 

1 

Silicon  Oil 

1 

10 

Glycerine 

1-25 

1000 

Preliminary  to  various  tests  with  liquids,  runs  were  rsaae  with  empty 
gyroscope  in  order  to  establish  the  nutations!  damping  rate  due  to 
bearing  friction.  The  damping  was  always  smal3,  about  .017  per  second, 
ana  was  applied  as  a  correction  In  all  tests. 

A  characteristic  of  the  tests  with  liquids  was  that  instability  was 
self-generating  only  in  the  vicinity  of  resonance.  The  width  of  the 
resonance  band  within  which  self-generation  occurs  (either  in  terms  of 
±  percent  of  fill,  or s )  increased  with  higher  spin,  higher  specific 
gravity  of  the  liquid,  and  lower  j  values.  The  effect  of  spin  can  be 
seen  in  Figure  3-  For  c/a  =  3*09  end  100  percent  fill,  the  instability 
is  seif -generating  at  6000  rpm;  at  4000  rpm  an  artificial  initial 
disturbance  was  necessary. 

For  the  c/a  =5-16  cavity,  the  combined  experimental  results  and 
theoretical  prediction  are  shown  in  Figure  4.  Figure  5  shows  the  results 
for  oil  and  glycerine  at  4000  ana  6CC0  rpm  in  greater  detail.  Similarly, 
Figures  6  and  7  show  the  experimental  results  for  c/a  =  3-09- 

As  can  he  seen  from  Figures  4  through  7-  the  resonance  band  at 
shorter  frequencies,  i.e.,  higher  fill-ratios,  is  cut-off  at  about  tq  = 
0-30,  producing  a  truncated  distribution.  To  test  whether  the  resonance 
band  is  symmetrical  about  its  maximum,  a  slightly  shorter  cavity,  c/a  = 
3.03,  was  tested.  Shortening  the  cavity  has  the  effect  of  shifting  the 
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FIG.  6 
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resonance  to  lower  fill -ratios.  The  results  are  shown  in  Figure  8.  For 
glycerine,  the  band  is  still  truncated,  but  for  oil,  the  response  is 
symmetrical  about  its  maximum. 


4.  DISCUSSION 


Stewartson's  theory  of  the  instability  of  liquid-filled  shell  indicates 
that  in  the  vicinity  of  coincidence  (resonance)  between  one  of  the  fluid 
frequencies,  xQ,  and  the  nutational  frequency  of  the  shell,  t  ,  the 
amplitude  of  the  nutational  component  of  yaw  will  grow  as 

Ojt 

e  (or^  =  cut  ) 

where  w  is  axial  spin  of  the  shell  and 


,  *  1/2  J  S  -  (t0  -  t/ 


s  = 


pjggl2^ 

I  a  c/a 
x  ' 


"Stewartson's  parameter" 


p  =  fluid  density 

28  is  a  tabulated  function  to  be  found  in  Stewartson's  tables; 

it  varies  with  x  and  percent  of  fill, 
a  =  radius  of  the  cavity 

o=7  1-1/s,  s  =  gyroscopic  stability  factor. 


For  instability,  x  must  be  real,  or 

f 


(x. 


>  0 


which  leads  to  the  well  known  Stewarts on  instability  criterion: 

uii ■  win*  'Wo 


(to  ~Tn>  , 

57s  < 
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GYROSCOPIC  SPIN  RESULTS  -  SHORTENED  CAVITY 


The  half-width  of  the  hand  of  instability  is 


At  =  /S 

The  maximum  rate  of  divergence  is  at  exact  resonance,  tq  =  t^  ,  for 
which 

(a,  )  =  ^  /S  per  second. 

v  1'max  2  r 

The  theoretical  curves  were  computed  using  these  relations. 

Observed  results  both  in  free  flight  and  with  gyroscope  differ 
significantly  from  the  predicted  behavior.  The  differences  appear  to  be: 

(a)  observed  maxima  of  divergence  are  considerably  lower  than 
predicted  and  show  marked  dependence  on  vis cosily 

(b)  the  width  of  the  resonance  band  also  is  considerably 
broader  than  predicted;  the  width  increases  with  viscosity 

(c)  there  appears  to  be  a  shift  of  (t  )  relative  to  t 

‘  v  o'max  n 

with  viscosity. 

It  is  instructive,  as  a  first  guess,  to  correlate  these  effects  with 
Reynolds  numbers. 

4.1  Correlation  with  Reynolds  lumbers: 

We  define  the  Reynolds  number  of  the  fluid  in  the  cavity  as: 


The  following  are  the  Reynolds  numbers  of  various  tests: 


v  c.s.  Re 

? 

Free  Flight,  Glycerine  1000  ?.0  x  10 

Gyroscope: 


4,000  rpm 

6,000  rpm 

Mercury 

0.1 

- 

6.3  x  106 

Water 

1 

4.2  x  105 

6.3  X  lo5 

Oil 

10 

4.2  x  10* 

6.3  X  10* 

Glycerine 

1000 

4.2  x  1C2 

6.3  x  102 

We  shall  consider  the  three  observed  differences  in  turn. 
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a.  The  Height  at  Resonance 

The  ratio  of  theoretical  to  observed  rate  of  divergence  at 

resonance 


x  obs 


is  plotted  vs-  Re  in  Figure  9*  Gyroscope  results  correlate  quite  veil. 

From  the  trend  of  the  correlation  curve  it  would  appear  that  an  agreement 
with  theoretical  prediction  would  be  attained  at  higher  Reynolds  numbers. 

An  attempt,  therefore,  was  made  to  test  mercury.  Tests  with  c/a  =  3-09 
ana  3-03  cavities  were  unsuccessful  because  of  violent  divergence.  Partial 
success  was  achieved  in  the  c/a  =  5-16  cavity.  The  results  are  shown  in 
Figure  10.  An  extrapolation  of  the  results  to  resonance,  an  uncertain 
pr  jcedure,  suggests  that  the  resonance  would  occur  near*  t^,  as  predicted, 
and  the  rate  of  divergence,  at  resonance,  would  be  close  to  the  theoret¬ 
ically  predicted  value.  However,  if  one  incorporates  into  Figure  9  two 
other  available  data  points:  frcro  free  flight  tests  and  Ward’s  result 
(see  Appendix  M,  these  disagree  with  gyroscope  results. 

b.  The  Width  of  the  Resonance  Band 

The  resonance  band  shows  progressive  widening  as  viscosity  in¬ 
creases.  At  higher  viscosity  it  also  exhibts  surprisingly  wide,  almost 
non -terminating,  tails.  However,  along  these  tails,  the  instability  of 
the  gyroscope  is  not  self -generating  but  requires  a  bit  of  artificial 
help  to  start  the  nutational  motion  going.  Once  disturbed,  the  motion 
is  divergent .  With  an  empty  cavity  a  similar  initial  disturbance 
invariably  damped. 

With  such  broad  tails  it  is  impractical  to  define  the  width  of  the 

band.  For  our  purpose  we  define  the  width  at  one  half  of  its  maximum 

amplitude  and  compare  this  with  theoretical  width  at  corresponding 

height,  i.e.,  v  =wi  .  at  l/2  of  (or.)  .  This  ratio,  so  defined, 

—obs  J-th  —  max 

is  plotted  vs.  Re  in  Figure  11.  Again,  there  is  good  correlation  of 
gyroscope  results  but  disagreement  with  free  flight  and  Ward’s  values . 
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c.  Displacement  of  with.  Viscosity 

Examination  of  the  resonance  graphs  shows  clearly  that  there 
is  a  progressive  displacement  of  the  observed  peak  of  the  resonance  curve 
with  viscosity.  Thus,  with  glycerine,  the  peak  occurs  at  tq  =  .026 
(c/a  =  3-03)  whereas  the  nutational  frequency,  =  .040,  remains  unchanged. 
In  other  fineness  ratios.  5.36  and  3-C9>  the  peaks  are  similarly  displaced. 
However,  in  these  cases  the  positions  of  the  peaks  for  glycerine  are  quite 
uncertain  and  are  emitted  from  the  plot.  The  peak  of  other  liquids  also 
might  he  somewhat  influenced  by  considerable  asymmetry  of  the  resonance 
curves . 

The  observed  difference,  -  Tn  is  plotted  vs.  Re  in  Figure  12. 

In  this  case,  the  formerly  discordant  two  data  points,  one  from  free 
flight  and  another  from  Ward,  are  in  somewhat  better  agreement  with 
gyroscope  results. 

The  question  may  well  be  raised  whether  as  modified  by  the 
inertial  properties  of  the  liquid  (Appendix  3)  is  the  appropriate  refer¬ 
ence,  or  should  not  xn  of  the  empty  system  be  used  instead?  With 

*o 

mercury,  for  example,  the  observed  t^  =  .043  whereas  that  of  the  empty 

gyroscope  is  Tn  =  .036.  The  resonance  peak  appears  much  closer  to  .043 
.  o 

than  to  .03o.  Also  for  glycerine  the  peak  is  at  ?Q  =  .026,  whereas  the 

observed  t  =  .040  and  that  of  the  empty  gyroscope  is  tR  =  .038.  The 
n  'o 

displacement  is  clearly  present.  For  the  case  of  free  flight  firings, 

the  observed  peak  is  at  t  =  .055?  whereas  t  is  .061  and  with  liquid 

o  nQ 

.065.  Thus,  the  displacement  of  (t  )  appears  to  be  real. 

o  max 

4.2  Viscous  Damping: 

Relatively  good  correlation  of  the  gyroscope  results,  in  the  form  of 
r,  w,  and  At  with  Reynolds  numbers  indicate  the  Importance  of  viscosity. 
However,  a  significant  disagreement  between  these  results  and  r,  w,  and  Ar 
as  obtained  from  the  free  flight  and  from  Ward's  experiments  suggest  that 
simple  correlation  with  Reynolds  numbers  is  inadequate. 


27 


To  examine  the  effect  of  viscosity  through  the  Navier-Stokes  equations 
is  very  difficult-  However,  in  general,  the  effect  of  viscosity  on  the 
oscillating  fluid  system  is  to  act  as  an  energy  sink  or  a  damper.  In  a 
mechanical  analog  this  effect  would  he  representable  by  a  dash-pot.  The 
solution  for  damped  fluid  oscillations,  therefore,  would  contain  imaginary 
components  in  all  resonant  frequencies  terms.  Any  fluid  frequency,  there ^ 
fore,  would  be  of  the  form 


-r  =  t  +  iS 
o  oo 


(I) 


where  ?  is  the  resonant  freauency  of  an  inviscid  fluid  and  6  is  the 
oo 

damping  factor  of  the  tq  mode. 

From  reference  3  cr  5>  the  nutational  root  of  the  characteristic 
equation,  defining  the  oscillations  of  the  system,  can  be  written 


T 


(2) 


where  S  is  the  "Stewartson  parameter”  previously  defined.  The  solution 
of  this  guadratic  for  t  leads  to 


i  = 


i  +  T 
o 


-O2  • 


(3) 


The  nutational  amplitude,  in  the  absence  of  aerodynamic  damping,  is  given 
by 

imTt 

constant  x  e 

For  instability  only  the  imaginary  part  of  t  is  important.  Moreover,  we 
must  have 

3  -  (To  -  \)Z  >  0  • 

which  is  Stewartson’ s  instability  criterion,  and  take  the  negative  sign 
in  front  of  the  square  root. 

With  viscous  damping.  Equation  (2)  becomes 


T 


T 

n 


/i.\ 
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The  solution  for  x  becomes: 


T  "  *■  I  [5  *  4'  &2  ‘  (Too  -  Tn)2  ‘  1  26<Too 


Let 


m=S+52-(x  -  t  )2 

oo  ny 


n  =  25(x  -  t  )  . 

oo  ny 


The  real  negative  part  of  the  imaginary  component  of  complex  t  is: 

r~w 

/  m  *  ym  + 


i M 


-  5 


(6) 


At  resonance,  i.e.,  when  xqq  -  x^  =  0,  Equation  (6)  becomes  particularly 


simple.  Call  this  viscous  x 


xv  =  1/2 j^  +  62  -  bj 


(7) 


The  inviscid  x.  is 

i 


x.  =  1/2  Js 


a.  The  Amplitudes 


(8) 


The  problem  now  is  to  determine  5  from  the  experiments.  We 
assume  that  the  observed  ratio  of  amplitudes,  r,  is  equal  to  the  ratio  of 
inviscid  to  viscous  x*s,  i.e.. 


r  = 


th  Ti 


JS_ 


°1  obs 
From  Equation  (9) 


Tv  /s  +  b2  - 


(9) 


5  =  1/2  JS  (r  -  i) 


(10) 


(5) 
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The  following  table  summarizes  the  values  of  r,  w,  and  At  previously- 
plotted  in  Figures  9>  H,  and  12.  In  the  ease  of  At  only  the  averages  are 
given. 


Observed  r,  v,  At 


Source 

/s 

r* 

V 

-At  x  10? 

Range 

1.37  x  10"2 

3.7 

2.3 

10 

Ward 

0.87  x  10'2 

4.4 

3.3 

5 

Gyroscope: 

4000  rpn 

6000  rpa 

Hg 

13.69  x  10‘3 

- 

(1.4) 

- 

(1) 

(0) 

K?0 

3.71  x  10~3 

2.8 

- 

1.4 

- 

2 

Oil 

3-71  x  10*3 

5.9 

5-6 

5.0 

2.9 

5 

Glycerine 

4.16  x  10-3 

15.1 

11.9 

8.7 

8.1 

16 

For  the  gyroscope,  the  values  of  fs  are  those  for  c/a  =  3.O9. 

x 

A  plot  of  log  6  vs.  log  Re  suggested 
_  const. 

O  =  ■  — "  • 

(Re)n 

3ecause  of  the  distribution  of  the  data,  neither  the  constant  nor  the 
value  of  n  could  be  well  determined.  However,  the  following  values 
represent  the  data  reasonably  well. 

n  =  l/4 

constant  =  0.135  ±  -017  s.d. 

A  comparison  of  computed  t./t^  ("predicted")  with  the  observed  r  is  shown 
in  Figure  13a.  There  are  some  indications  that  at  lower  Reynolds  numbers, 
for  Re  <  lO5,  the  exponent  n  night  be  l/p  and  be  more  nearly  l/2  at  higher 
Reynolds  numbers.  3ut  present  data  are  insufficient  to  make  this  distinction. 

The  form  of  dependence  of  §  on  Reynolds  numbers  vas  suggested  by 
Weaemeyer  on  Theoretical  grounds ,  see  Reference  6. 
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FIG.  13 


b.  The  Width 


In  contrast  to  the  inviscid  solution  which  exhibits  definite 
band  width,  2/s,  the  viscous  solution.  Equation  (6),  shows  tails  whose 
amplitudes  only  asymptotically  approach  zero*  Hence,  it  is  impractical 
to  define  the  resonance  band  width.  However,  a  very  simple  expression 
for  characteristic  width  can  be  obtained  if  we  neglect  m2  relative  to 
n  under  the  square  root  of  Equation  (6).  Such  characteristic  width, 
then,  can  be  defined  by  the  following  equation! 

-6=0.  (11) 

From  this  it  follows  that  the  viscous  half -width,  is 

wy  =  /S  +  S  .  (12) 

The  inviscid  half -width  is 


The  ratio  w  =  w ^/w  =  1  +  b/JSy  or 

6  =  /S  (w  -  1)  .  (13) 

The  width,  w^.  Equation  (12)  occurs  from  one  third  to  one  h»lf  of 
the  maximum  amplitude.  This  depends  on  the  value  of  6//S.  Therefore, 
the  ratio  w  =  w^/w^  is  not  strictly  comparable  to  the  previously  used 
ratio  w  at  one  half  of  the  maximum  amplitude.  Nevertheless,  the  two 
are  compared  in  Figure  15b.  For  "predicted"  values,  w ^/w^,  previously 
found  6  =  -155/  v'Se  was  used.  The  agreement  is  fair. 

c.  Displacement  of  (t  ) 

office 

Interaction  between  oscillating  fluid  particles  and  the  walls  of 
the  cavity  apparently  gives  rise  to  the  boundary  layer.  This  boundary 
layer,  in  turn,  may  produce  an  "effective"  fineness  ratio  or  the  cavity, 
slightly  different  from  the  geometrical.  A  rough  estimate  of  this 
effect  can  be  made  as  follows. 


Let  c/a  be  the  geometrical  fineness  ratio.  If  the  boundary  layer 
has  thickness  y,  the  "effective"  ratio  will  be 

g..  -  y 

a  -  y 


The  difference  is 


A(— )  =  s.  .  2-zjl 

a  a  a  -  y 


\a  a 

1-1 

a 


From  Stewartson’s  tables  one  finds  that  a  change  A(— )  produces  a  change 

a 

At  .  The  two  axe  related  by 
o 

A(f)  =  1.2(2j  +  1)  Atq  .  (15) 

Equating  Equations  (l4)  and  (15),  the  desired  answer  is 


At  = 
o 


(1-1) 


Since  At^  is  known  from  observations,  (1)  can  be  computed  and  plotted  vs. 
Re-  The  results  are: 

A  .  c/a  -  1  .086 

Ar  =  T  -  T  - - iy:  ■-  X  -  •  (17) 

o  oo  n  2j  +  1 

A  plot  of  "predicted"  Atq,  computed  by  Equation  (17),  vs.  observed  is 
shown  in  Figure  15c.  The  agreement  is  tolerable. 

Finally,  it  is  of  interest  to  compare  the  observed  resonance  band 
with  that  theoretically  predicted  with  viscous  damping  included.  The 
theoretical  curve  was  computed  using  Equation  (6),  i.e.. 


(o,)M  -S 

-  theory  n\ 


/  /  2  2 
u)|  /  m  +  Jm  +  n 


Vr5 
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and  compared  with  the  observed  cn’s  for  oil.  Figure  8,  c/a  -  3*03,  6000 
rpm.  The  comparison  is  given  in  Figure  14.  The  agreement  is  excellent. 
Particularly  noteworthy  is  the  theoretical  appearance  of  wings  in  the 
resonance  band  which  were  so  conspicuous  in  gyroscope  experiments  with 
more  viscous  fluids.  Therefore,  introduction  of  viscous  damping  clearly 
improves  the  inviscid  theory  when  applied  to  real  fluids. 


5.  CONCLUDING  REMARKS 


Stewartson*s  inviscid  theory  still  provides  the  best  available 
guide  for  predicting  a  priori  the  dynamic  behavior  of  liquid  carrying 
shell  with  cylindrical  cavity.  It  is  suggested,  however,  that  shell 
designers  will  do  better  by  taking  into  account  the  effects  of  viscosity 
of  their  fluids.  The  following  steps  outline  the  procedure  to  be 
followed. 


a.  ilutational  frequency:  This  to  be  computed  by  taking  inertial 
properties  of  the  liquid  into  account. 


I 

Tn  =  i/2  -£  (i  +  a) 

y 

where 


I 


x 


+ 


(18) 


(19) 


I 

y 


I„  ,  I  are  the  axifl  and  transverse  moments  of 

--o  yo 

the  empty  shell  and  ix  and  iy  are  those  of  the 
rigid  body. 


inertia  respectively  of 
fluid  regarded  as  a 


b.  The  maximum  divergence  will  occur  at  the  principal  fluid  frequency 


given  by 


where 


Jmax 


t  -  At 
n 


Ar 


c/a  -  1  .086 

2J  +  1  Jsfii 


(20) 


(21) 
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c.  From  Stewart son’s  tables  one  finds  the  value  of  (c/a)/(2j  +  l) 

corresponding  to  above  ~  .  Let  this  be  K  .  Therefore, 

°max  o 

(c/a)Q  =  KQ(2j  +  1)  . 

These  are  the  fineness  ratios  of  the  cavities  at  which  max-iwurn  instability 
will  occur. 

d.  Because  of  the  width  of  the  resonance  band,  the  designer  should 

avoid  fineness  ratios  in  the  vicinity  of  the  critical  (c/a)  .  Half -width 

o 

of  the  resonance  band  is 


max 


=  /S  + 


5 


Therefore,  one  should  find  from  the  tables  the  Undoing  fineness  ratios 
corresponding  to 

To  =  To  *  C/S  +  8)  *  (22) 

°  max 

Roughly,  this  is  given  by 

A(|)  =  1.2(2  j  +  l)(/§  +  6)  .  (2?) 

Therefore,  the  fineness  ratios  tc  be  avoided  are 

(|)  ±A(|)  .  (24) 

o 

It  should  be  noted  that  the  inviscid  instability  criteria 


-i  <  <  i 


/ s 


with  viscous  correction  becomes 


wnere 


(25) 
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r 


=**,i 


An  example:  Prom  our  range  firings,  =  .065.  From  Equation  (21), 

At  =  .011.  Therefore,  the  maximum  divergence  will  occur  at  fluid  frequency 

T°mar  =  ~  -011  =  *054.  From  Stewartson*s  tables  for  b2/a2  =  0.10 

(90  percent  fill)  we  find 


zfTi  =  1.044  =  Kq 


(c/a)o  =  5.15 


For  this  shell, 

Js  =  1.57  x  10"2 

8  =  2.65  x  10"2 
JS  +  8  =  .040 

and  by  Equation  (24),  the  range  of  fineness  ratios  to  be  avoided 
5.13  ±  .14 

or  from  c/a  =  2-99  bo  c/a  =  5.2?.  From  Figure  2  and  tables  we  obtain 
the  following: 


Lower  limit 


Unoer  limit 


t  =  .020 

o 

To=  .090 


c/a  =  5.09 
c/a  =  5.27 


These  numbers  are  in  a  fair  agreement  with  the  estimates  by  using 
Equation  (25)- 
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With  slightly  greater  sophistication,  the  width  of  the  resonance 
band,  or  the  range  of  (c/a)  values  about  the  critical  should  be  determined 
relative  to  aerodynamic  damping.  Only  the  rate  of  divergence  of  the 
nutational  component  of  yaw  which  is  in  excess  of  the  rate  due  to  aero¬ 
dynamic  damping  is  important  in  practice.  If,  therefore,  a  is  the 
aerodynamic  damping  rate,  then  the  width  of  the  resonance  band  is  to  be 
found  from  the  condition 


(26) 


From  this  one  finds  lower  and  upper  and,  from  the  tables,  corresponding 
fineness  ratios. 
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APPENDIX  1 


PHYSICAL  AND  DYNAMIC  CHARACTERISTICS  OF  20mm  MODELS 

-  models  were  fired  in  this  program.  Since  their  inferred  characteris- 
hoved  no  systematic  variation  with  fineness  ratios  of  the  cavities, 
were  combined.  All  errors  assigned  to  the  mean  values  are 
rd  errors. 

oral  Frequency  cp-  '  =  .0271  £.0001  rad/cal 

ricr.al  Frequency  o2  '  =  .0060  £.0001  rad/cal 


-Y 


o.  '  +  a,  '  =  .0331  ±.000h  rad/cal 
v  =  .-25  rad/cal 

4s-  =  .'77 '''7  £.0011 


G  =  ^L—  —  ?  ,  =  .61-0  £.012 

3,  SL 


*  It  <lw> 


.065  £.001 


rom  physical  measurements  of  empty  models  prior  to  firings. 


with  the  above  z,  rives  t_  =  .0-1.  The  difference  between  eirmtv 

o  *  * 

r.d  tine  dynamic  4s-.  as  inferred  from  firings,  is  due  to  the  effect  of 
•-rtial  properties  of  the  liquid. 


Also: 


ana 


C  =  3.66  ±  .06 

O' 


1  I  s  v2 

S  =  ~-^T  =  1.69  ±  .02  =  rT-xf 


i-o2 


2IYp  S  d3C, 


M 

or 


Kncwinr  the  ra*es,  one  can  detennine  Ix  and  IY  separately. 
Le* 

I.  =  I.  +  o'  i.  _ 


IY  =  IY()  +  3  iyo 


where  ix  and  iY  refer  to  the  liquid  acting  as  a  rigid  body.  We  have: 


IX{?  =  2li.8l  gm-crri2  i„o  =  2.7^  gm-cnr2 

ly c  =  332.09  gn-cm2  iyQ  =  17.68  gm-cm2 

The  firings  indicate  that 

a  =  .96 

3  =  -70 


All  firings  v:sre  done  at  2000  fps  from  a  gun  with  a  twist  of  rifling 
1:15.  Other  aerodynamic  characteristics  of  the  models  are: 

CD  .350  Wt  gms 

C.  2.17 

if 

a 


C;  '  -.012 


k2 


Damping  ra'es: 


Uutation 

ttj  -  3-20  x  10'4 

per  caliber; 

F recession 

=  2-5l:  10“' 4 

eer  caliber 

-5 


APPENDIX  ? 


THEORY  OF  GYROSCOPE 

This  theory  is  sufficiently  well  known  not  to  require  full  exposition 
ere?  see,  for  example,  reference  7.  However,  a  few  brief  remarks  may  be 


Let  I,  be  the  polar  moment  of  inertia  of  the  rotor 

Iv ,  I_  I..."  be  the  transverse  moments  of  inertia  of  the  rotor, 

inner  ana  outer  gimbals  respectively,  about  the  y  axis 
I.  ,  I.  '  similarly  about  2  axis.  Clearly,  Iz "  =0. 

Let  lyo  =  Iy  +  V  +  V' 

I;c  =  Iz  +  Iz' 

if  the  gyres cooe  is  so  designed  that 


-  c 


X.  leu, 

I  =  M..  I-  be  the  moment  about  either  axis. 
i  c  -  o 

Let  ur  assume  that  bearing  friction  can  be  represented  as  f79y  and 
3-  respectively.  Moreover,  if  f,  =  f.  ,  as  was  the  case  with  our 
bearinrs,  then  the  differential  equations  of  motion  can  be  written  as: 

I  9..  +  f  Sy  +  I,  n0.  -  K09y  =  0 

I  3-  -i-f  S;  -  lx  nsy  -  K0 9,  =  0 

•where  n  is  constant  spin  in  rad/sec,  and  Kq  =  gl-Ih  is  gravitational  moment. 

M  is  the  mass  of  the  rotor  and  the  inner  gimbal  and  h  is  the  center  of  mass- 
pivot  separation.  For  the  unstable  position,  i.e.,  when  c.m.  is  above  the 
nivot .  ?L,  >0. 


(1) 


b5 


f 

1 


Define:  9 


0Y  +  i  02 


The  two  equations  can  be  combined  to  give 


where 


I  9  +  (f  -  i  Ixn)9  -  Ko0  =  0 

i  =  A  . 


"he  characteristic  equation  of  (2)  is 


I  m2  +  ff  -  i  Ixn]m  -  K0  =  0 


The  roots  are  complex  giving  two  frequencies  and  two  damping  rates. 
With  an  approximation  that  f/lxn  «  1,  these  are 


Nutation 


.frequencies 
Ixn  .  . 

IT 


Damping  Rates 

-  k 


Pi*0C0Ssion 


I*n  r-,  i 

? r  (1-o) 


(l-o  ) 


c  -  4 1 '  - 


s  =  r=? —  gyroscooic  stability  factor. 
Non-dimensional  freauencies  are 


n  21 


(l-*o) 


%  - 1  ^ 


If  o  ~  1,  s  is  large,  Tn  =  =4-  and 


o=!-^+  ... 


_  I  .  Ko 

P  21  2s  Ixn5 


It  should  be  noted  thau  if  Kc  <  0,  statically  stable  position,  s  <  0, 


a  >  1  and  the  preeessional  amplitude  has  a  tendency  to  diverge  at  a  rate 
depending  on  the  quality  cf  the  bearing  or  on  the  value  of  f.  For 
statically  unstable  configuration ,  K0  >  0,  s  >  0  and  o  <  1,  friction 
damps  both  amplitudes. 

If  w‘i  introduce  into  one  component,  say  9» ,  an  additional  stabilizing 
torque,  -K,  due  to  measuring  device,  the  equations  become 

I  9v  +  Ixn  02  -  (K0  -  K)9V  =  0 
I  9,  -  I,n  9,  -  K0  e;  =0 

where  friction,  -which  does  not  affect  the  frequencies,  was  neglected. 

It  can  be  shown  that,  in  this  case,  the  new  frequencies  are: 


'  n 


I* 

I 


as  before,  and 


APPENDIX  3 


CONTRIBUTION  OF  LIQUID  TO  THE  INERTIAL  PROPERTIES  OF  THE  SYSTEM 


Both  in  free  flight  and  gyroscope  experiments  with  liquids  the 
observed  nutational  frequency  was  always  somewhat  higher  than  for  the  empty 
system.  The  liquid,  therefore,  must  have  contributed  its  inertial  properties 
to  the  system.  From  free  flight  experiments  it  was  found  that  the  liquid 
contributed  96  percent  of  its  axial  moment  and  ~0  percent  of  its  transverse 
moments  of  inertia  when  considered  as  an  equivalent  rigid  body. 

The  gyroscope  offered  an  opportunity  to  find  the  liquid’s  contribution 
with  greater  accuracy.  The  obvious  choice  of  fluid  is  mercury  for  which 


Again, 


we  define: 


I  =  I0  +  3  i0 

where  IXq,  Tj,  and  ix  ^ ,  i0  are  the  axial  and  transverse  moments  of  inertia 
of  the  rotor  and  of  the  liquid  considered  as  a  rigid  body,  respectively. 
From  Appendix  2,  it  is  easily  seen  that 


^9 


K 

I 


(Tn 


+  T  V* 
P' 


n 

Si 


and 


G 


7n-TP 

Tn+Tp 


The  freauencies  are  directly  observable  Quantities.  If  one  knows  I.  , 

o 

and  Iq  and  the  moment  K  from  prior  physical  measurements,  then,  by 
experiments,  one  determines  Ix  and  I,  and,  hence,  a  and  3.  A  series  of 
experiments  were  run  for  this  purpose  and  the  results  are  given  in  the 
followinr  table: 


Dynamic  Properties  of  Mercury-Filled  Gyroscope 
c/a  =  5*16 
Spin  I*. 


%  Fill 

rpm  x  10  3 

I 

or 

3 

100 

2.5 

.OU37 

.n 

-69 

4.0 

.OUio 

-98 

.82 

90 

4.0 

.0454 

•  97 

.74 

65 

3-0 

.0456 

-97 

.67 

4.0 

.0451 

•  95 

.68 

6.0 

.0446 

.9^ 

.69 

So 

3-0 

.0456 

.06 

.63 

75 

3-0 

.0458 

.95 

.67 

Average 

.96 

.70 

zrinp'CY 

.0359 

> 

Ithough  the  exact  agreement  of  average 

a  and  3  with 

free  flight 

value . 

is  fortuitous,  the 

experiments  indicate  that  the  a; 

■rial  moment 

9O 


of  the  empty  system  is  fully  augmented  by  the  axial  moment  of  the  liquid, 
considered  as  a  rigid  body,  and  the  transverse  moment  by  TO  percent. 

The  latter  value  appears  to  be  in  good  agreement  with  the  theoretical 

Q 

value  ior  an  inviscid  fluid  .  In  the  inviscid  theory,  there  is  no 
contribution  by  the  liquid  to  the  axial  moment  of  inertia  of  the  system. 
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APPENDIX  1* 


WARD’S  EXPERIMENTS 


Ward  5  tested  Suewartson's  theory  with  a  gyrostat  by  varying  fill- 
ratios  in  a  fi::ed  cylindrical  cavity. 

It  is  to  be  recalled  that  Stewartson’s  instability  condition  is: 

(To_Tn) 

-1  <  - t -  <  1 

S2 

Ward  found  that  the  following  limits  best  satisfy  his  observations: 
o  „  _  (To"Tn^ 

-3-9  < - t —  <  2.7 

S2 


Thus,  the  limits  are  considerably  broader  than  predicted  and  asymmetric 
relative  to  r  .  However,  he  did  find  that  maximum  instability  appeared 
to  coincide  with  the  predicted  position  at  66  percent  fill-ratio. 

Frca  the  cited  reference,  the  characteristics  of  the  gyrostat  are: 
Cavity:  a  =  I.U3  cm 

c/a  =  3. 00 
Volume  =  55  cc 
Hutationai  frequency,  Tn  =  .112 
Fluid:  v  =  23.9  c.s. 

p  =  1.2  gm/cc  (assumed) 

Spin  6,000  run 
Re  =  5.U  :<  Id3 

'With  above  characteristics  of  the  cavity  and  =  .112  ve  find,  from 
Stewartson’s  tables  at  t.  =  t_  ,  2R  =  .30S  (Stevartson * s  R).  Therefore 


55 


I 

s”2  =  11L.3 

Thus, 

-3.9  <  ll"o(70_7=  )  <  2.7 
it  follows: 


70 

b3  /a2 

£  fill 

Lower  limit 

.07  8 

.28 

72 

Upper  limit 

.136 

• 

LO 

CD 

62 

Average 

-107 

•33 

67 

According  to  these  estimates,  if  the  resonance  band  is  sjrmmetrieal  about 
its  maximum,  the  maximum  should  occur  at  t0  =  .107  or  at  67  percent  fili- 
ra:io.  This  irrolies  a  shift  of  maxima.  from  i 

n 

.107  -  -112  =  -.005 

Such  a  displacement  of  maxima  -will  remove  the  asymmetry.  In  our  experiments, 
see  Figure  8  for  oil,  the  resonance  curve  is  symmetrical. 

The  width  of  the  band  At  =  .058,  whereas  theoretical  width  =  .017. 
The  ratio,  therefore,  is  3.3.  Moreover,  free  Ward’s  graph  of  the  rate  of 
divergence  of  the  nutational  amplitude  at  resonance,  the  amplitude  doubles 
in  about  1.1  seconds.  In  our  notation,  this  implies  <r2  =  0.6.  The 
theoretical  value  Is  jjJS  =  2.7b.  The  ratio,  therefore,  is  h.u. 

These  results  were  used  in  comparison  with  our  experiments . 


5^ 


